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ABSTRACT 

X- In this paper for the quark Nambu-Jona-Lasinio model and for the 

qj \ phenomenological four fermionic model of QCD we analyze vacuum energy 

^ ■ structure and calculate dependence of energy of one-particle excitations 

on momentum. We reveal existence of two different types of excitations 
above the vacuum with the respective minimization of the vacuum energy 
rS | density. First type of excitation energy spectrum is the usual spectrum of 

j^ ■ relativistic massive particle, but the second one has a linear dependence 

on momentum. 
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1. Introduction 

In their pioneer paper 1 Nambu and Jona-Lasinio (NJL), using the analogy with su- 
perconducive mechanism, had shown an intimate connection of unitary-nonequivalent 
transformations of quantum states with the existence of nontrivial vacuum. Later 
this approach was generalized on a case of quark interaction in chiral QCD limit 
(see, e.g., review 2 ). In this paper we consider one particle excitations of vacuum and 
its energy spectrums for wide set of permissible transformations. We show that be- 
sides the known NJL solution there always exists one more excitation with a linear 
spectrum on momentum and this fact is a property of any contact four -fermionic 
SU(2) x SU{2) x t/(l)-symmetric interaction. 

2. Analysis of the NJL solution. 

The NJL Hamiltonian is 



^ d r 

H = J d 3 xip(x)p^fip(x) / d 3 x 



2 /-, , Cj . , ,\2 



(ip(x)if>(x)\ — (if>(x)^f 5 rif>(x 



where p — —i\J, color and flavor indexes on if) are assumed. The NJL transformation 
representing Heisenberg field if>(x) at t — in terms of massless spinors u(p), v(p) 
leads to the following form of the Heisenberg field (dynamical mapping 3 ) 



^(x,0) = j / d 3 p N r (p)a r (p)e ipx + M r (p)b +r (p)e- ipx (2) 

N r (p) = (cose{p)+ 1 ° sin 6{p))u r {p) 
M r {p) = (cose(p)- 1 °sme(p))v r {p), 

u r (p), v r (p) are the massless spinors ($u(p) = ftv(p) = 0) normalized by the usual 
conditions u +r (p)u s (p) = v +r (p)v s (p) = 5 rs . It is a simple exercise to check that the 
spinors N r (p) and M r (p) satisfy the equations: 

(E(p)f -p^-m*(p))N r (p) = 

(E(p)-f -prf + m*(p))M r (p) = 0, (3) 

with E(p) = p/ cos26(p), m*(p) = pt&n26(p), p—\ p |. Creation and annihilation 
operators a r (p), b r (p) and a +r (p), b +r (p) satisfy the conventional anticommutation 
relations for fermions and are interpreted in a sense of dynamical mappings as oper- 
ators of excitations of "physical" particles. Vacuum state is defined with respect to 
them. 



Vacuum energy Wo — (0 | H | 0) evaluated by using Eqs.(|^) is turned out to be a 
functional of the transformation angles 9(p): 



V 3 (2vr) 



- N r N f G 



c^/l 



(27T) 



: , , kcos(26(k))d 3 k- 



sm(26(k))d 3 k 



2G ± N C N 



I 



(V*) 



(4) 



where G = G x {\ + 2N c N f ), ± = j^js J d 3 k, V = J d 3 x. Thus, 9(p) "enumerates" 
all possible vacuums. Variation of 8(p) gives the following condition for the extremum: 



G 



(2k) 



d 6 k 



m 



pta.ia29(p) = const. 
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m* ^0 



(5) 



Another trivial solution corresponds to m* = 0. To answer the question whether 
the solution suits for maximum or for minimum of Wq we need to check the sign of 
functional determinant: 



det 
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56(k)56(p) 
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J F 
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d 3 k 
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E(kf 
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G 



(27T) 



d 6 k 



m 



E(kf 



> 0. (6) 



N 



^Jd 3 kln(8N c N f E(k)) 



So, if m* ^ this extremum is minimum which is degenerated at m* = into the 
inflection point. 

The energy spectrum of one-particle excitation E*(k) is determined by the follow- 
ing eigenvalue equation: 



H,d +r (k)} | 0) = E*(k)d +r (k) |0) 



(7) 



d +r (k) stands for a +r (k) or b +r (k). Calculating the commutator in the l.h.s. of above 
equation using the dynamical mapping (0) and omitting three-particle excitations we 
obtain expression for E*(k): 

E*(k) = k cos(26(k)) + — sin 2 (26(k)). (8) 

On the solution (|[) E*(k) coincides with E(k) 

E*{k) = Vk 2 + m* 2 = E{k). (9) 

Now we want to draw attention to the important fact, that under the functional 
variation of Wo over 6(k) one more solution is lost, and it should be studied separately. 
Namely, the solution: 8(k) = = const. This solution, as will be shown below, admits 



an extremum in Wq, and this extremum corresponds to minimum of Wo- In this case 
m* is a linear function of momentum: m*(k) = ak, and the expression for the vacuum 
energy has the form: 



(10) 
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Jkd 3 k 

<k> = fffik- 


3m ([II]) we obtain: 


dcj) V 


<fc> , / G 


cos 20^ = 


oil! Ziiy i 

k* V < fc > v* 


1/v - i««yC*>»o. 


.(|TT|) determines the parameter in the extremum point: 


^ <k>V* , 
COS 20 = — = c , 





(11) 

(12) 

(13) 

and this extremum is in fact always minimum (Eq.(^)). Hereafter to the end of the 
text we will deal with positive values of c*, so that the minimum existence condition 
is just restriction c* < 1 (recall, that we use the system c = h = 1) and this means 
c* is less than speed of light. Inserting the expression of cos 20 via c* (Eq. fll3|) ) at 
6{k) = into Eq.(||) we finally obtain the energy spectrum E* of excitation (0): 

E*(k) = c*k + ^(l-c* 2 ), (14) 



As one can recognize from Eq.(14), c* has the meaning of group velocity of the 
excitation. Thus, the spectrum is proved to be a linear function of momentum , 
i.e. "neutrino" type, and the condition of the minimum existence corresponds to the 
demand c* < 1. At c* = 1 the minimum of Wq is degenerated into the inflection 
point. 

So, NJL transformations for the Hamiltonian (P yield to two types of stable ex- 
citations above vacuum. One excitation describes a massive particle with the normal 
relativistic spectrum, and the mass arising dynamically. The second excitation has 
the linear energy spectrum. This fact is a consequence of unitary-nonequivalenceness 
of the transformations (0) from the massless particles to the massive ones causing 
the vacuum to be nontrivial. As will be shown later appearing of these two types 
of excitations is not a result of the given transformation (0) and of the Hamiltonian 
choice ([!]). They arise generally for any contact four-fermionic SU(2) x SU(2) x U(l)- 
globally symmetric interaction. Let us note, eventually, that the NJL transformations 
can be represented in terms of quantum fields as: 



^(x,0) =e G g(x,0)e 



-G 



G= j d 3 xd 3 y[q(-\x,0)£(x- y)g (_) (x,0) - g (+) ( x , 0)^(x - y)g (+) (x, 0)] , (15) 

where q(x) = g < - + - ) (x) + g < -~- ) (a:) is massless spinor field expanded on its frequency parts, 
and 6(k) is Fourier image of the real function £(x). 



3. Model of pseudo-scalar couplings. 

Consider now SU(2) x SU(2) x £/(l)-invariant Lagrangian 2 : 



ipiM+ % 



ip^i^Tip^x] 



2 
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ip(x)ip(x] 



iI)(x)t~i^(x)) + (ip(x)-f 7 M r0(a;) 



ip{x)a lll/ ip{x)\ + [if>(x)ia liv 'y ti/j(x 



(16) 



Our aim is to quantize the respective Hamiltonian by the dynamical mapping method 3 
in such a way that the relevant vacuum would contain condensate with the quantum 
numbers of 7r-mesons. Therefore, we take generator of the transformations from 
"physical" fields to the Heisenberg ones (dynamical mapping) exploiting the expres- 
sion for the axial charge. 

^(x,0) = e%(x,0)e- 6 (17) 

where q p (x, 0) is a "trial" physical massive field with a mass m(k) arbitrarily depend- 
ing on momentum. 

g«(x,0) = -^T,jd 3 kg(k)^-u r (k)a r (k)e i ^, 



(2tt 



g(~)(x,0) = *£ / ' dPkg(k)^v*(k)b»(k)e 



ikx 



u r (k) and v r (k) are massive spinors normalized by the conditions u r (k)u s (k) 
—v r (k)v s (k) = 5 rs , E{k) = wk 2 + m?{k). Operators a r and b r anticommute as: 



{b r (k), b +s (p)} = {a r (k), a +s (p)} 



E{k) 
m(k) 



S rs S(k-p)I c I f , 



(19) 



I c and If are unit flavor and color matrixes (tr/ a = N a ). Wave function g(k) is 
introduced for " smooth" regularization of divergent integrals and to attach physical 
meaning to them. If quantums of physical fields are interpreted as in paper 3 , then 
g(k) characterizes a momentum distribution inside the excitation, and the integral 
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(2^W / I d(k) I 2 d 3 k = V* determines V* - space volume of this excitation. We 
normalize | g(k) | 2 in such a manner that it is ~ 1 in the region where it essentially 
differs from zero, i.e. we choose: 



9(k) 



1 1 g(k) | 4 d 3 k 
1 1 g{k) | 2 (Pk 



(20) 



Correspondence with the usual momentum "cut-off" parameter A is established due 

to the equality: 

I j-A i 

So, in the local (over momentum) relations we will put | g(k) | 2 = f but saving it in 
the integrals to regulate its convergence. Taking into account all above comments we 
construct the generator G from the diagonal (with respect to the frequencies) parts 
of the axial charge: 



G 



d 3 x 



,Av5^R 



g^( x , OhVr^"^, 0) - g w (x, 0) 7 u 7 5 fg w (x, 0) 



a(+)< 



^5- (+), 



(21) 



where f = £ ■ (nr), £ = const, n is a unit vector in the flavor space. With such 
choice of G we obtain the following dynamical mapping: 



^(0,x) 

N r (p) 
M r (p) 



\Y.!d'k m{k) 



27r)iV^ E{k) 

cos(£§(k)) + inT>y >y 5 sm(£$(k)) u r (p) 

cos(M>(/c)) - mr7°7 5 sin(£$(/c))| v r (p) 



g{k)N r {k)a r {k)e lkx + g{k)M r {k)b +r {k)e 



-ikx 



(22) 



From above it follows that £ and $(fc) = ^4 play the role of parameters of G- 
transformation, i.e. as in Eq.(^) they "enumerate" all possible unitary- nonequivalent 
vacuums. For the vacuum energy density we find out the expression: 



(0 | H | 0) 
V 



G 

~2 
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d 3 k | g(k) 



E{k) 



cos{2£${k)) + 



n 2 



d 3 k\g(k) | 2 ■$(*;) cos (2£$( A;)) 



(2vr) 3 
-^ f d 3 k\g(k) | 2 ■sin(2£$(A;)) 



+ 



(27T) 



+ const, 



(23) 



where G = G\ + 2G\N c Nf + 2G3, and the const does not depend on £ and $(&). 
Hence it follows, firstly, that the terms of the Lagrangian (|T6| ) corresponding to the 
vector and axial-vector hadronic currents (p and A\ mesons) do not contribute to 
the vacuum energy density. Secondly, each SU(2) x SU(2) x [/(l)-invariant term of 



the Lagrangian ([16]) gives the uniform contribution to Wq, so that accounting all the 
terms of ([TB]) results in the summary constant G. Let us vary Wq on <&(k), reckoning 
the physical vacuum to be a minimum of the energy density W /V, and equating the 
variation to zero. Using the relation 6(k. 2 /E(k)) = —m(k)5§(k) we have an integral 
nonlinear equation for m(k) 

1 k + Aii • *{k)) tan(2l$(A;)) = & + ^ ~™ {k) (24) 



\E(k) ^ V 'J v K " ^ 2£ 

/ii = J^Lj d3k I *K*0 I 2 "$(A:)cos(2£$(A;)); 

»2 = j^Jd 3 k\g(k)\ 2 - S m(2£$(k)). 

This equation should be true for the all admitted values of k: < k < oo. The 
analysis shows that the equation (24[) has the only solution m(A;) = afc, with a = const 
satisfying the transcendental equation: 

1 / 2£a \ a 

tan -===--. (25) 



VT+a 2 \VT+a?J 2£ 

This equation is thought to connect the parameter £ with the constant a. There are 
many solutions of fl25|) written as 



2£a / 3\ / 1 , 

2n(n + - ! + — + ••• , n = 0,l,2,... (26) 



VTT^ V 47 V a 2 

The expression for Wq/V under condition ( p5| ) has the following form: 

f = -2^1^ ( /(0) _ ^Z—fta)) + »n*, (27) 

denoting 

cos ( 7 f^) 

/(a) = 7==r , ( 28 ) 

VI + a 

and 

/fc|fl(fc)| 2 rf 3 fc 

9 / I 9(k) I 2 d 3 * • 



Now, in order to fix the value of £ (or a) we have to minimize fl27|) over £(a). The 
very form of ([27]) hints to take derivative with respect to a rather than £, so we do: 

MV - *«W^ (A.))' >0 (29) 



d W 


d 2 W 



Thus, as in above, the extremum exists and always corresponds to minimum. First 
of Eqs.(|D links f(a) with c* 

f(a) = ^<k> g =c\ (30) 

besides the expression for c* via the parameters of the model at hand coincides exactly 
with the expression obtained under the NJL transformations Eq.(|H|), while the second 
of Eqs.(P§D constrains c* < 1. 

Evaluation of excitation energy E* from Eq.(0) results (compare with (|j)) 

k 2 

E* = — — cos(2£0(k)) + ui<$>(k) cos(2£$(k)) + a 2 sm(2£$(k)) (31) 

E{k) 



What is most surprising, the Eq.(pl|) exactly coincides with expression fll4|) on the 
solution fl30"D 

E*(k)=c*k + ^(l-c* 2 ). (32) 

If in Eq.(^) the parameter £ is chosen to be a function of k : £ = £(k), that 
corresponds to the input of a function £(x — y) into the generator G Eg. (|2T|) , then 
it is not hard to make sure that the equation on the extremum of vacuum energy in 
this case yields for m(k) the NJL solution, i.e. m(k) = m* = const. 

We had studied all possible generators G, composed from the bilinear combinations 
of fields. From the results of calculations there follows the only conclusion that in the 
models with the contact four fermionic interaction possessing the SU(2) x SU(2) x 
U(l) symmetry two types of stable excitations emerge - with the universal expressions 
for one particle energies Eq.(^) or Eq.(p!4|) and with respective (in terms of m* or c*) 
universal expressions for vacuum energies. The linear spectrum at c* = 1 corresponds 
to the neutrino spectrum. However, this point cannot be achieved because of two 
reasons. 

First: at c* = 1 the minimum degenerates into the inflection point. 

Second: rewriting (on m(k) = ak with 1/a = go) the expression for the spinor 
u r (k) in the form 

1 




u r (k) = , , =■ l + Vl + g 2 -go7 K ( 33 ) 



(k = nk, q = ng ) 
and finding out the following behavior of go near c* ~ 1 



?o ~ n ^ (34) 



Vi 



c 



*2 



one can see, owing to the relativistic factor in Eq.([34]), the second obstacle to achieve 
the point c* = 1. 



Finally note, that the relativity - invariant generators G, constructed from all 
bilinear combinations of fields can be brought to any of the forms a + b + , ab, ect. How- 
ever, when the algebra of the generators is observed it is infinite-dimensional because 
of the relativistic phase volume d 3 km(k)/E(k), since every generator contains differ- 
ent orders of m(k)/E(k). The exception concerns the case of linear spectrum with 
m(k)/E(k) = const, then algebra is closed and the parameter space becomes finite- 
dimensional. 

It is worth mentioning in a conclusion that the energy spectrum can be rewritten 
in the more convenient form for the experimental treatment, namely: 

E*(k)= X ock + -^(l-xl). (35) 

Here Xo is a dimensionless number, c is the speed of light in "emptiness", and the 
vacuum stability condition reduces in this case to Xo < 1- Parameter y^ has dimension 
of energy and, as well as the parameter xo, is not determined by the initial Lagrangian. 
Its determination needs a model for the distribution function g(k). In general these 
parameters should be extracted from experiments. If E* is identified with spectrum 
of physical neutrino, then we are to suggest the following (probably not complete) 
list of existing neutrino experiments, where this spectrum has chance to be found: 

a) Neutron beta-decay (spectrum of electrons near kinematical bound) 

b) Tritium beta-decay (Tr — > He 3 + e + i? e ) at the end of electron energy spectrum, 
where at present there exists discrepancy between theory and experiment. 

c) Double neutrinoless beta-decay (test on the Majorano mass). 

d) Solar and atmospherical neutrino (oscillations in "vacuum" and matter). 

e) Neutrino signal from the supernova (time delay and oscillations) 

Authors are deeply grateful to Prs. D.A.Naumov, S.S.Gerstein, V.A.Matveev and 
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